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ABSTRACT 

Thermal conduction in tangled magnetic fields is reduced because heat conducting 
electrons must travel along the field lines longer distances between hot and cold regions 
of space than if there were no fields. We consider the case when the tangled magnetic 
field has a weak homogeneous component. We examine two simple models for tem- 
perature in clusters of galaxies: a time-independent model and a time-dependent one. 
We find that the actual value of the effective thermal conductivity in tangled magnetic 
fields depends on how it is defined for a particular astrophysical problem. Our final con- 
clusion is that the heat conduction never totally suppressed but is usually important in 
the central regions of galaxy clusters, and therefore, it should not be neglected. 

Subject headings: magnetic fields: conduction — cooling flows — magnetic fields: dif- 
fusion — methods: analytical 



1. Introduction 

In order to solve the problem of electron thermal conduction in a stochastic magnetic field, 
one should consider separately two effects that reduce the conduction (Pistinner & Shaviv 1996; 
Chandran & Cowley 1998). The first effect is that the heat conducting electrons have to travel 
along tangled magnetic field lines, and therefore, they have to travel longer distances between hot 
and cold regions of space (Tribble 1989, Tao 1995). The second effect is that electrons, while they 
are traveling along the field lines, become trapped and detrapped between magnetic mirrors, regions 
of strong magnetic field (Chandran, Cowley, & Ivanushkina 1999). 

In our recent paper (Malyshkin & Kulsrud 2000) we studied the second effect, and we found the 
reduction of the thermal conductivity parallel to the magnetic field lines, Ky, relative to the Spitzer 
value for the thermal conductivity. Kg, caused by the presence of magnetic mirrors. In this paper 
we consider the first effect, and we calculate the further reduction of thermal conduction relative 
to At||, caused by the tangled structure of the magnetic field lines. As a result, we obtain the total 
reduction caused by the both effects, and we calculate the effective thermal conductivity Keg, which 
one has to apply for astrophysical problems. It turns out that there is no an universal definition of 
KcfT) and the result depends on the particular astrophysical problem under consideration. 

The structure of this paper is the following. In the next section we consider the random 
walk of tangled magnetic field lines by making use of a diffusion approximation for this random 
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walk. We find the expressions for the probabihty distributions of the field line lengths and for 
the Laplace transform, in the field line length, of these distributions. In Section 3 we use the 
Laplace transforms to calculate the effective thermal conductivity for our stationary model of a 
galaxy cluster. In Section 4 we consider a time-dependent model and study the time evolution of 
temperature in a galaxy cluster. Finally, in Section 5 we discuss our results, compare them with 
the results reported in previous papers, and consider the applications to cooling flow problem. For 
two opposite limiting cases our results for Kqs are: one is similar to that of Tribble (1989), and the 
other coincides to that of Tao (1995). 



2. Diffusion approximation for the random walk of magnetic field lines 

In a particularly simple approach, the behavior of tangled magnetic field lines can be considered 
by using the one-dimensional random walk model suggested by Tribble (1989). In this model the 
lines are assumed to random walk between two infinite boundary plates, which are placed at x = 
and X = Xq perpendicular to the a;-direction. In this section we consider this simple model in 
order to find the probability distributions of the lengths of the field lines in clusters of galaxies 
(we then use these distributions to find the effective thermal conductivity in the next sections). 
Contrary to the discrete calculations of Tribble, we use the continuum diffusion approximation for 
the random walk of the field lines. This allows us to include a weak homogeneous mean magnetic 
field component into our calculations. 



2.1. Diffusion equation 

We assume that the mean magnetic field component, (B), is homogeneous, i.e. (B) = const. 
Let choose the coordinate system in a such way that (B) lies in the x-z-plane: {B)^ = {B) cos 9, 
{B)y = 0, (-8)2 = (-B)sin^. Here 9 is the angle between the mean field component and the x- 
direction. Further, assume that the random component of the magnetic field, 5B, is much stronger 
than the mean component, i.e. {B)/SB <C 1. Let express x-, y- and z-components of 5B in the 
spherical system of coordinates: 6Bx = SBcos^p, 5By = 5B simp cos and 6Bz = SB sin ip sin ip. 
The random component is isotropically distributed, so cos (p and ip are uniformly distributed over 
[—1,1] and [0, 27r) respectively. Under these assumptions, the cosine of the angle between the total 
magnetic field, B = 5B + (B), and the x-direction is 

cos a = cos (p + [{B) I bB\ [cos 9 — cos 9 cos^(p — sin 9 sin Lp cos 95 sin V'] , (1) 

where we keep terms only up to first order in {B)/6B <C 1. 

Let the largest scale of the random component of the magnetic field be Iq <^ Xq. The largest 
scale component has more magnetic energy in it than all smaller scale components have. Therefore, 
the decorrelation length of the total random component is Iq, and over each decorrelation length 
the random field changes into an entirely new direction. As a result, the mean step and the 
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mean-squared step of the field line random walk in the x-direction are 

(Ax) = loicosa)^^ = {2/3)lo€Cos9, 



(2) 



over each decorrelation length. Here, to obtain these final expressions, we average cos a and cos^a 
over direction and over absolute value of the field random component (i.e. over cp, and SB). We 

def 

also introduce the parameter e = {B)/{SB) <^ 1, which is the mean ratio of the magnetic field 
mean component to the random component. 

There is the unique magnetic field line that goes through any given point of space between 
the two boundary plates. This field line leaves the point along two branches: a positive branch 
that starts at the initial point and goes always in the direction of the local field, and a negative 
branch that goes always opposite to the direction of the local field (see Figure 1). Let consider the 
following problem. Start at point (x = X, y, z) and follow the positive branch of the magnetic field 
line going through this point, i.e. along the field line and always in the direction of the local field. 
The positive field line branch random walks in space according to equations (2). Let P^{s, x) ds dx 
be the probability that we are at x-position x G [x,x + dx) [and at any y-, z-positions] after we 
"have walked" along the line a distance s G [s, s + ds). The upper index "(+)" refers to the positive 
branch of the magnetic field line. Using equations (2), it is straightforward to write the diffusion 
equation for P~^{s,x) as ^ 

dP+ _ {Ax)dP+ 1 ((Ax)2) _ 2ecosedP+ lod^P+ 

ds Iq dx 2 Iq dx^ 3 dx 6 dx'^ ' 

P+(0,x) = 5{x-X). (4) 

This is a Fokker-Planck equation. Here, S{x — X) is the Dirac delta-function, and the formula on 
the second line is the initial condition on P~^{s,x), which means we start walking along the field 
line at X = X. If we now consider the field line's negative branch, which goes always opposite to 
the direction of local field, then we need to replace e by — e in formula (3) in order to obtain the 
differential equation for the probability function P~{s, x) for this negative branch. (Note, that this 
replacement corresponds to the substitution vr — 6* for 6*.) 

Hereafter, we use convenient dimensionless variables 

s={SX^/lo)~\, x = x/Xo, X = X/Xo, P± = XoP±. (5) 
Note that the positions of the boundary plates are x = and x = 1. In these variables 

^-^"^-^^^ P-(0,x)=^(x-X). (6) 

Here, we introduce an important parameter 

/? =^2XoZo^ecos0. (7) 



^Note that the diffusion approximation is valid if a; S> /o and Xo — x ^ lo- 



-4- 



In the limit /3 ^ 1, the effect of the mean magnetic field component is negligible and the 
"diffusion term" (1/2) ((?^P^/5x^) in equation (6) is dominant. In this small /3 limit the random 
walk of field lines is controlled by their diffusion in space, and for a given s ~ 1, the functions 
P''=(s,x) have considerable spreads in x. On the other hand, in the limit /? » 1, the effect of 
the mean magnetic field component is large and the "flux term" —(3{dP^/dx) in equation (6) is 
dominant. In this large f3 limit the probability functions stay very narrow in x for a given value 
of s, P^(s,x) « d(x — X ^ /3s). When /? ~ 1, both terms in the right-hand side of equation (6) are 
equally important. Thus, /? indicates whether the influence of the mean field component on the 
field line random walk is negligible (/? ^ 1), important (/3 ~ 1) or dominant (/? S> 1). 



2.2. Boundary conditions and probability fluxes 

Equations (6) need the boundary conditions on the probability functions P^(s,x) at the two 
boundary plates. Let assume for definiteness, that the plate at x = represents the boundary 
of a galaxy cluster, the plate at x = 1 represents the cluster core (see Figure 1), and the mean 
component of the magnetic field has a positive x-component, i.e. < ^ < 7r/2 (clearly the sign 
of cos 9 has no bearing on thermal diffusion problem). As we said above, there are two field line 
branches that we can follow starting at x = X, the positive branch and the negative branch. Each 
branch random walks in space, and can reach the cluster boundary at x = and the cluster core 
at X = 1. 

Let first consider the boundary conditions at the cluster core. The core is assumed to be much 
denser and much colder than the rest of the cluster, so we can neglect the temperature T in the 
core and set it to zero. As a result, whenever the positive or the negative branch reaches the core 
at X = 1, it looses its thermal energy and cools down to zero temperature immediately. Therefore, 
there is no any point in following the field line branches after they first reach the core, and the 
corresponding boundary conditions on the probability functions P^ at the core are the absorption 
conditions: P'''(s, 1)=0. 

Let now consider the boundary conditions at the cluster boundary. We assume that the 
density drops significantly there. Therefore, the random component of the magnetic field also drops 
considerably at the cluster boundary, because this random component is believed to be created by 
MHD dynamo action inside the cluster. On the other hand, the mean component of the field hardly 
changes at the boundary, because it is associated with the open magnetic field lines that leave the 
cluster. As a result, the parameter /3, given by equation (7), increases and becomes large at the 
cluster boundary (while the mean field component can still be less than the random component, so 
that our diffusion approximation is still valid). In other words, we consider (3 to be constant inside 
the cluster and to become large at the boundary. Therefore, at the cluster boundary x = 0, the 
term —(3 (dP^/dx) in equation (6) is always dominant. As a result, the positive field line branch is 
reflected at the cluster boundary, while the negative branch is absorbed and leaves the cluster. 

To summarize, we write the boundary conditions for the positive and negative field line 
branches as 

P~(s, 0) = 0, absorption of the negative branch at the boundary, , . 

P~(s, 1) = 0, absorption of the negative branch at the core; 
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/3P'^(s,0) — (l/2)(5P"'"/5x)|j^^Q = 0, reflection of the positive branch at the boundary, , . 

P+(s, 1) = 0, absorption of the positive branch at the core, 

see Figure 1. 

Let for a moment consider the positive field line branch and integrate the diffusion equation (6) 
over X G (0, 1), using boundary conditions (9). As a result, we obtain the negative change of the 
total probability that we are still inside the cluster, x G (0, 1), after we walked a distance s along 
the line's positive branch: 



P+ s,x dx 



ds ./o ' 2 dx 



= F^(s). (10) 

=1 



Here we introduce a "probability flux" into the cluster core at x = 1: 

F+(s) = -(l/2)(5P+/5x)|,^, (11) 

(the upper index "+" refers to the fleld line's positive branch). Note, that the probability flux 
across the cluster boundary at x = is zero because of the reflection condition there [see eqs. (9)]. 
The probability, that we leave the cluster and enter the cluster core when s G [s, s + ds), is equal 
to (s) ds. Therefore, (s) is the probability distribution of the lengths s of all positive field line 
branches that start at x = X, random walk in space and finally come to the cluster core at x = 1. 

Now consider the negative field line branch. In this case, we integrate the appropriate diffusion 
equation (6) over x G (0, 1) and use the boundary conditions (8). We find that the negative change 
of the total probability is now given by the sum of two probability fiuxes F (s) and Fq (s) into the 
core (at x = 1) and across the boundary (at x = 0) respectively: 

Fr(s) = -(1/2) idP-/dx)\^^„ Fo-(s) = (1/2) {dP-/dx)\^^,. (12) 

Similarly to the case of the positive branch, Fj^(s) and Fq (s) are the probability distributions of 
the lengths of all negative field line branches that start at x = X, random walk in space and finally 
come to either the cluster boundary at x = 0, or to the cluster core at x = 1. 

Because the probability fluxes are the probability distributions of the lengths of field line 
branches, hereafter, we refer to them as to the probability distributions. 



2.3. Laplace transform solutions 

Although equations (6) are linear, they are still difficult to solve analytically for both sets 
of boundary conditions given by equations (8) and (9). We solve equations (6) numerically in 
Section 4. However, in some cases we do not need the solution of equations (6) in order to find 
the effective thermal conductivity in tangled magnetic fields. In the next section we will find this 
conductivity using a simple stationary model and the Laplace images of the probability distributions 
Fj'"(s), F]^(s) and Fq (s) that we obtain in this section. 

To calculate these Laplace images, we first take the Laplace transforms s — > s, P^(s,x) —>■ 
?=*"(§, x) of equations (6). We have 

/9P± 1 «2p± 

sP±-i(x-X) = + (13) 
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The Laplace images P^(s,x) must be continuous functions of x. Thus, integrating equations (13) 
across x = X, i.e. over an infinitesimal interval x G (X — 0,X + 0), we obtain an additional jump 
conditions together with the continuity conditions 



Here, we denote the jump at x = X as [...]x±o- 



[P 



±1 



Jx±o 



0. 



(14) 



The boundary conditions on the images P (s, x) are the same as those on the original functions 
P^(s, x). In the case of the negative branch they are given by equations (8), in the case of the positive 
branch they are given by equations (9). Equation (13) is a simple linear homogeneous differential 
equation provided x 7^ X. There are two independent solutions of this equation expressed in 
terms of hyperbolic functions. We express the two general solutions of this equation in the regions 
< x < X and X < x < 1 as two linear combinations of these hyperbolic solutions. Then, we 
use the appropriate boundary conditions and jump conditions (14) to find the unique continuous 
solutions of equations (13) for P=''(s, x) in the whole interval of x, < x < 1. We have 



P+(s,x) 



2 exp [I3{x - X)] 



[Co cosh (^oX) + /3sinh (^qX)] sinh [^0(1 - x)], X < x; 



^0 ^ocosh^o + /3sinh,^o [ [Co cosh (^qx) + /3 sinh (^ox)] sinh [^o(l - X)], x<X 
for the positive field line branch, and 



P-(S,X) 



2 exp [-/3(x - X)] 



Co sinh Co 
for the negative branch. Here, we introduce 



sinh (CoX) sinh [Co(l-x)], X < x; 
sinh (Cox) sinh [Co(l -X)], x<X 



(15) 



(16) 



(17) 



Now, we substitute these formulas into Laplace transformed equations (11) and (12) to find 
the Laplace images F]^(s), F^(s) and Fq (s) of the probability distributions F^(s), F^(s) and Fq (s). 
We have 



Fr(s) 

Fr(s) 

Fo-(s) = / 
^0 







e-^= F+(s)ds 



e-^=Ff (s)(is 



19P^ 



2 dx 
IdP- 



2 9x 
Id?- 



00 

^"'^Fo(s)cis = 



x=l 



x=l 



x=0 



Co cosh (CoX) + 3 siuh (CqX) ^;3(i_x) ^g^ 
Co cosh Co + /3 sinh Co 



sinh (CqX) g_^(i_x) 
sinh Co ' 
sinh[Co(l-X)] ^gx_ 
sinh Co 



(19) 
(20) 



Here, the integrals are the definition of the Laplace images, and Co is given by formula (17). 



As is well known from the theory of the Laplace transform [and can easily be checked by 
differentiation of equations (18)-(20) with respect to the Laplace variable s and setting s to zero], 
the following formulas stand for the integral moments of the probability distributions (s) , F]" (s) 
and Fq (s): 



Jo 



s)ds = (-l) 



n '1 



s=0 



/ s"F-(s) 
Jo 



ds 



s=0 



(21) 
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For example, the zeroth integral moment of F^(s) is 

F+(s)ds = F+(0) = 1, (22) 



/o 

which means that the total probability for the positive branch to ultimately reach the cluster core 
is unity, as it should be because positive branches always end up in the core. The zeroth integral 
moments of F]"(s) and ^q{s) are 



oo oo 



/ 





which means that F^ (0) and Fq (0) arc the probabilities for the negative branch to reach the cluster 
core (at x = 1) and the cluster boundary (at x = 0) respectively. Note, that the total probability 
for the negative branch to reach either the core or the boundary, is equal to one, F|f(0) + Fq (0) = 1, 
as it should be. 

It follows from equations (23) that Fj"(s)/F^(0) and F^(s)/F^(0) are the normalized (to one) 
probability distributions of the lengths of only those negative field line branches that reach the 
cluster core and only those negative branches that reach the cluster boundary respectively (see 
Section 2.2). The probability distributions of the positive branches, F^(s), does not need to be 
normalized, because all positive branches always end up in the cluster core. 



3. Stationary model 

In this section we consider a stationary one-dimensional temperature distribution in a cluster 
of galaxies. We assume that there is a stationary homogeneous source of heat q between the two 
Tribble boundary plates, which represent the cluster core and the cluster boundary. The heat is 
transported by electrons along the tangled magnetic field lines. We further assume that the thermal 
conductivity parallel to the field lines, is constant. We keep our assumption made in Section 2.2 
that the cluster core is cold, i.e. T = at x = 1. On the other hand, the density drops significantly 
at the cluster boundary, so we use the heat reflection condition there, dT/ds = at x = (here, s 
is the distance coordinate along a field line) . 

Let consider a point (x = X, y, z) inside the cluster. The stationary temperature at this point 
depends on the lengths of the positive and the negative branches of the magnetic field line passing 
through this point. The positive branch always reaches the cluster core at x = 1, where it is cooled 
to zero temperature. Let the positive branch have dimensionless length s^. As for the negative 
branch, there are two possibilities. First, with probability F^(0) it reaches the core at x = 1 and 
cools down to r = 0. In this case we denote the dimensionless length of the negative branch by 
sj~. The second possibility is that the negative branch with probability Fq (0) reaches the cluster 
boundary at x = 0, where dT/ds = 0. In this case we denote the dimensionless length of the 
negative branch by Sq • In any case, the differential equation for the temperature distribution along 
the field line is 



+ -2^-^ = 0, (24) 
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where s is the dimensionless coordinate along the field line, given by equations (5). 



Now, let us consider a plane x = X, which is perpendicular to the x-direction, and let us find 
the temperature averaged over points of this plane. We make this averaging by randomly choosing 
points of the plane, and then, by averaging temperature over the chosen points. ^ The lengths s^*", 
and Sq of positive and negative branches, going through different points are uncorrelated, and 
these lengths have the probability distributions (s) , (s) and Fq (s) respectively. Therefore, the 
averaged temperature and the averaged temperature squared are 



(T) 



9X 



2kii 



^ {Fr(0)(s+)(sr) + Fo (0) [{{stf) + 2(s+)(so-)] } , 



9X1 



(25) 



2kii 



/2 



{Fr(0)((s+)2)((sr)') + Fo (0) [((s+)^) + 4((s+)3)(so-) +4((s+)2)((so-)2)]}26) 



Here, the factors multiplying by F]~(0) represent the averaged temperature obtained by solving 
equation (24) in the case when both the positive and the negative branches reach the cluster core 
(where T = 0), as shown by the lower field line in Figure 1. The factors multiplying by Fq (0) 
represent the averaged temperature obtained by solving equation (24) when the positive branch 
reaches the cluster core and the negative branch reaches the cluster boundary (where dT/ds = 0), 
as shown by the upper field line in Figure 1. Remember that Fj~(0) and Fq (0) are the probabilities 
for the negative branch to reach the core and the boundary respectively, while the positive branch 
always reaches the core. All brackets (...) mean averaging over appropriate probability distributions 
of lengths of branches. Note, that the expressions inside the brackets {...} in the equations above 
are dimensionless, while q, /cy, Xq and lo are not. 

Now we use equations (21)-(23) to express averaged powers of the branch lengths in the 
formulas (25) and (26) in terms of the Laplace images of the probability distributions. We obtain 



(T) 



(T2 



9Xq1 



2kii 



I'o 1 
9X^ 



+ 



ds ds 



+ 2 



dFp dF+ 
ds ds 



5=0 



2kii 



/2 



d^Ft 



+ 4- 



.dF; 



d^F+ d2p-j2p+ ^^2p-^2p+ 



ds dH^ dp' ds^ 



(27) 



(28) 



s=0 



We can substitute the Laplace images of the three probability distributions given by formulas (18)— 
(20) into these equations and obtain analytical formulas for (T) and (T^). However, the resulting 
expressions are too complicated to be usefully interpreted. Therefore, we calculated all derivatives 
in the equations above numerically. Figure 2 shows the relative dispersion of the temperature, 
((T^) — (T)2)/(T)2, as a function of parameter (3 for two choices of x-position inside the cluster, 
X = 0.5^0 and X = O.IXq. We see that the dispersion is high for small values of /3, while it is 
oc <C 1 when /3 3> 1, according to the discussion in the last paragraph of Section 2.1. 



^Our description for the temperature distribution along the field lines really refers to a space average of temperature 

over a length interval long compared to the field decorrelation length, lo- The effect of regions of weak and strong 
magnetic field (magnetic mirrors) is taken into account by introducing fty, which is reduced relative to the Spitzer 
conductivity, see Malyshkin and Kulsrud 2000. 

^Formally, this plane is infinite, so are the number of points and distances between the points. 
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Fig. 1. — There is the unique magnetic field hne that goes through any given point (x = X, y, z) of 
space, this field line leaves the point along the positive branch and along the negative branch. All 
positive branches arc reflected at the cluster boundary at x = 0, and arc absorbed at the cluster 
core at X = 1. All negative branches are absorbed both at the boundary and at the core. Thus, 
there are two types of magnetic field lines: first, there are lines that go from the cluster boundary 
to the cluster core (see the upper line in this figure), and second, there are lines that leave the core 
and finally come back to the core (see the lower line in this figure). 
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Fig. 2. — The relative temperature dispersion, ((T^) — (r)^)/(r)^, versus (3. The solid and the 
dashed lines represent X = O.bXo and X = O.IXq cases respectively. 



-10- 



If there were no magnetic field at all, the temperature dispersion would be zero and the 
stationary temperature at x = X plane would be 

Ts = ^iXl-X^), (29) 

where the effective thermal conductivity would be equal to the Spitzcr thermal conductivity, 
Keff = us- Because there is tangled magnetic field, the actual stationary temperature is given 
by formula (27) and it is higher. However, instead of using formula (27) we can use the familiar 
formula (29) for the temperature if we choose the appropriate rcdiiccd vahic of the effective thermal 
conductivity. Equating formulas (27) and (29), we define the effective thermal conductivity for our 
stationary model as 

_ ll Xj-X^ j~ d^Ft , df^dft JfodftV .... 



s=0 



The dashed lines in Figures 3(a) and (b) show this effective conductivity normalized to k^^Iq/Xq as a 
function of /? for the two choices of position inside the cluster, X = O.IXq and X = O.SXq. To obtain 
the plots we again calculated all derivatives of the Laplace images in equation (30) numerically. 
Note, that is the parallel thermal conductivity reduced by random magnetic mirrors (Malyshkin 
and Kulsrud 2000). 

Note, that for large values of (3 the effective conductivity is independent of X and is simply 
Keff = '^||/3^^o/9"^o ~ '^11 (^^/^o)^- ^ This result exactly coincides with that given by the equation (4) 
of Tao (1995), who gives this result as a lower limit on KeS- On the other hand, for small and 
moderate values of f3 the effective conductivity depends on X as on a parameter. This dependence 
results from our definition of /Cgfr given above. In other words, the actual value of Keff depends 
on how one defines it for a particular problem under consideration. The calculations in the next 
section further support this important statement. 



4. Time evolution of temperature in clusters of galaxies 

In the previous section we considered a stationary one-dimensional temperature distribution in 
a cluster of galaxies, assuming a constant heat source inside the cluster. In this section we solve a 
time-dependent problem and find the evolution of temperature of the cluster in time. Let assume, 
that when the cluster was formed at zero time, t = 0, the temperature was homogeneous inside the 
cluster, i.e. T = Tq, 0<x< 1 (here x is the dimensionless x-coordinate inside the cluster; x = 
and X = 1 correspond to the boundary and the core of the cluster respectively, see Figure 1). We 



''This result can be found by the following calculations. The solutions of equations (6) in the limit /3 ^ 1 
are P*(s, x) w S{x — X =p /?s). Thus, all positive field line branches reach the core and have dimensionless lengths 

~ (1 — X)//3, while all negative branches reach the boundary and have lengths ~ X//9. Substituting these 
expressions and Fj"(0) ~ 0, F(j"(0) ~ 1 into formula (25), and equating the result with formula (29), we obtain 
Kes = k,\i(3^1q/9Xo = K\i{Ax/lo)^, see also eqs. (7) and (2). This result is valid only if P'^Iq/Xq <C 1, which is 
equivalent to the assumed condition e = (B)/(6B) <C 1, see Sec. 2.1. 
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assume that the cluster cools down in time by the heat conduction into the dense cluster core. As 
in the previous section, we again assume that the parallel thermal conductivity is constant, that 
the core is cold, T = at x = 1, and that there is the heat reflection condition along magnetic field 
lines at the cluster boundary, dT/ds = at x = (here s is the coordinate along field lines). 

Because electrons travel along magnetic field lines, each field line cools down in time individ- 
ually. Let us consider a temperature evolution at a point (x = X, y, z) inside the cluster. There is 

a single field line going through this point. Initially, at t = 0, the temperature at this point is Tq. 
Then, as the field line cools down, the temperature drops in time. It is convenient to introduce a 
dimensionless time variable 



/2 



t, 



(31) 



where p and Ch are the mass density and the heat capacity (per unit mass) of gas inside the 
cluster. For simplicity, we assume that the product pCn is constant. Then, the time evolution of 
the temperature at the point (X, y, z) is given by the following simple equations: 



dT{T, s) 1 d'^T 



dr 



9 52s' 



T(0,s)=ro, 



(32) 



where s is the dimensionless distance counted along the field line starting at the point (X, y, z), see 
equations (5). 

There are two possibilities for the field line going through the point {X,y,z). First, both the 
line's positive branch of length and the line's negative branch of length reach the cold cluster 
core (see the lower line in Figure 1). In this case the boundary conditions on the temperature 
distribution along the line are 

T(r, -sr) = 0, T(r,s+) = 0. 

With these boundary conditions the solution of equations (32) for the temperature at the point 
(X, y, z), where s = 0, is 



Tii(r,Sj;,s+) =ro- V(l + 2n)-iexp 

TT — i 



n=0 



9(sr+sJ 



+^2 



(1 + 2nf 



sm 



vrs 



lSi +Si 



^(l + 2n) 



(33) 



Here, the index "11" indicates that both branches reach the core at x = 1. The second possibility 
is that the line's negative branch of length Sq reaches the cluster boundary, while the line's positive 
branch of length s'l ends up in the cluster core (see the upper line in Figure 1). In this case we 
have 



idT/ds)\ 



= 0, T(r, s+) = 0. 



and the solution of equations (32) for the temperature at the point (X, y, z), where s = 0, is 



roi(r,So,s+) =To- V(l + 2n)-i 

TT — i 



exp 



n=0 



36(so +s+)2 



(1 + 2nf 



sm 



TTSV 



2(so +S+ 



-(l + 2n) 



• (34) 
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Here, the index "01" indicates that the negative branch goes to the boundary at x = 0, while the 
positive branch goes to the core at x = 1. 

Now, let find the temperature averaged over all points in the plane x = X in the same way as 
we did in the previous section. We have 

oo oo oo oo 

{TKr) = I y"rn(r,sr,s+)F-(sr)F+(s+)dsrcis+ + J J To,{t, ,s+)f^ {s^)Ft {3+) ds^dsf . (35) 











Here, we have averaged over the probability distributions of sf, and Sq , which are given by 
equations (11) and (12). The two integrals in formula (35) correspond to the two possibilities 
for the field line branches, considered above. We have carried out these integrals numerically. 
Figures 4(a) and (b) show contour plots of {T){t)/Tq in the r-/3 parameter space for the two 
choices oi X, X = 0.5Xo and X = Q.IXq. To obtain these plots, first, we numerically solved 
equations (6) with the boundary conditions (8) and (9), by making use of the standard implicit 
algorithm for partial differential equations. ^ Then we calculated the probability distributions (11) 
and (12) for different values of (3. Finally, we substituted these distributions and formulas (33), (34) 
into equation (35), and we calculated the temperature for different values of r and (3. (The Laplace 
transform method was not useful for this problem because we need more than simple moments of 
the si etc.) 

If there were no magnetic field, the time evolution of the temperature would be given by 
equations 



dT Keff d'^T 



T{0,x) = To, {dT/dx)U=o = 0, T{t, Xo) = 0, 



dt pCn d'^x 

resulting in the following expression for the temperature similar to equation (34): 



(36) 



Ts{t) = To- V(l + 2n)-iexp 

TT ^ — i 



n=0 
oo 



vr 



Kest 



To-V (l + 2n)-^exp 



n=0 



4 X^pCn 
4 11k. 



(1 + 2nf 
(1 + 2n)2 



sm 



sm 



TrjXo - X) 
2Xo 

Tr{Xo - X) 



2Xo 



(1 + 2n) 
(1 + 2n) 



(37) 



where wc use definition (31). The effective thermal conductivity KefT would be equal to the Spitzer 
thermal conductivity, = Kg, if there were no magnetic field. 

For each set of values T, r, /3 and X there is a value of Kcs, which makes equations (35) 
and (37) agree. This vahie can be taken as the effective thermal conductivity, but it does depend 



®We used a grid with equal steps As and Ax in s- and x-coordinates. The goal was to find the tabulated solution 
of equation (6) at all grid points. The solution at s = is a delta-function. To find the solution at s > 0, wc used a 
method similar to mathematical induction. This method allowed us, after we obtained the tabulated solution at s, 
to find the unknown tabulated solution at s + As by the following algorithm. First, we wrote equation (6) as a finite 
difference operators, taking all x derivatives at s + As (the implicit algorithm). As a result, we obtained a system 
of linear equations, which was a tridiagonal matrix. Second, we solved this system by the method of the Gaussian 
decomposition with backsubstitution, and found the solution at s + As. 
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Fig. 3. — The effective thermal conductivity Keff normahzcd to kjiZq/Xq is plotted as function of 
parameter /3. The function Keff(/3) is not universal and it depends on how Keff is defined. 
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Fig. 4. — Contour plots of the average temperature normalized to the initial temperature, 
{T){t)/Tq, for two choices of X: (a) X = 0.5Xo and (b) X = O.IXq. 
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on conditions. In the case of the presence of the tangled magnetic field, instead of using the 
complicated formula (35) we can use the simple formula (37) with this Kpfr reduced appropriately 
from the Spitzer value. To find the reduced effective conductivity, we first choose fixed values of 
position X and temperature T. Then, we find Kgff cLS cl function of /3 in a such way, that both 
formula (35) and formula (37) give the same chosen value of temperature, T, at the same time r. 
The calculation of HcsiP) was done numerically for two choices of X: X = 0.5Xo and X = O.IXq; 
and for three choices of T: T = 0.03 Tq, T = 0.50 Tq and T = 0.90 Tq. [In other words, we 
numerically found such functions Kee{l3), that the temperature contours T/Tq = 0.03, T/Tq = 0.50 
and T/Tq = 0.90 appeared the same on Figures 4 when we used formulas (35) and (37) for the 
temperature.] The resulting functions Kcfr(/5) normalized to k^^Iq/Xq are given by the solid lines in 
Figures 3. We again see, that Kes{P) is not an universal function, and that it varies for different 
choices of T and X, except when the /? is large. For smaller values of temperature the effective 
conductivity becomes less because very long field lines, which keep the initial temperature, become 
more important. 



5. Discussion 

We would like to start the discussion with by stressing one of the main results we found: 
the actual value of this effective conductivity in the tangled magnetic fields depends on how this 
effective conductivity is defined. In our paper we used two "natural" definitions of k^s for two 
very simple models of a cluster of galaxies: a time- independent model in Section 3 and a time- 
dependent model in Section 4. The results for KqS, reported in figures 3(a) and 3(b), are the 
same only provided /3 3> 1, and they are significantly different (by up to a factor of ten) when 
/? < 1. We conclude that there does not exist an universal result for the effective conductivity in 
the tangled magnetic fields. One has to define K^fi and to calculate it for a particular astrophysical 
problem that he/she considers. We believe that such calculation is possible for many problems, 
including numerical simulations of galaxy cluster formation and of cooling flows, by making use of 
our diffusion approximation method for the random walk of tangled magnetic field lines. 

It is useful to compare our results for the effective conductivity in the tangled magnetic fields 
with those reported in previous papers. As we said in the last paragraph of Section 3, in the limit 
(3^1 our result for coincides with the result obtained by Tao (1995). This simple formula (for 
P large) corresponds to lines being in order Xq/c in length, a result sometimes believed (remember, 
that e <C 1 is the ration of the field mean component to the field random component). But in 
fact, it is only valid as a lower limit on KeS, except for large (3. However, we can speculate that it 
revolves around the very definition of e. We give e first and then impose the statistics of the field. 
If the statistics is given first and the mean field is defined afterwards, as the rms value of the mean 
field on the scales ^ Xq, then Kgff would be closer to our results. 

In the limit /9 <C 1, our result, k^s ~ 0.1k^^Iq/Xq, is consistent with the result of Tribble 
(1989). ^ To the best of our knowledge, there have no results for k^s when /? ~ 1 obtained before. 



^Tribble has Kch — 4k[|/q/Xq, sec his equation (11). The diflerence arises because: first, our boundary condi- 
tions (8) and (9) for the random walk of field lines are different from those of Tribble; second, the magnetic field lines 
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Rccently Chandran and Cowley (1998) have suggested that the diffusion of the heat conduct- 
ing electrons perpendicular to the magnetic field lines is crucially important for the heat transport. 
In their model the statistically independent random step by the electrons is equal to \/2Lr^Io, 
where Lrr ~ loln{lo/pe) is the Rechester-Rosenbluth length. This length is the distance along 
magnetic field lines over which a separation between two initially neighboring lines grows expo- 
nentially from the electron gyro-radius pe up to the field decorrelation length Iq (because of the 
Kolmogorov-Lyapunov exponential divergence of field lines). Chandran and Cowley suggested 
that the perpendicular diffusion of electrons and this exponential divergence of the magnetic field 
lines enable the electrons to conduct heat between different field lines, which may enhance the 
effective thermal conductivity considerably. However, we believe that the picture considered by 
Chandran and Cowley is rigorously valid only if electrons can move a distance Lrr along field 
lines without collisions. In fact, there are collisions, and if the electron mean free path A is less 
than the Rechester-Rosenbluth length Lrr, the exponential divergence of field lines does not help 
electrons to diffuse perpendicular field lines. In clusters of galaxies, the ratio of Lrr to A is 
Lrr/X « 5000(/o/10kpc) (r/10^K)-2 (n/10-3 cm ). As a result, the perpendicular heat trans- 
port can generally be neglected (we accept 3-10 kpc as a typical scale for magnetic fields in clusters 
of galaxies, see Kronberg 1994, Eilek 1999). The Chandran and Cowley's model can be very im- 
portant for some other astrophysical problems, but the consideration of their model is beyond of 
the scope of this paper. 

Finally, let estimate the importance of the heat conduction for the formation and evolution of 
galaxy clusters, using figures 4. The parameter /? and the dimensional time r, given by equations (7) 
and (31), can be estimated as 



/3 = 20 



Xq 10 kpc ecosO 
IMpc ~To 0l~' 



r = 5x10 



-6 '^ll 



10 



Ks lO^'^yrs 



n 



5/2 



In 



10 kpc 



IMpc' 

Xn 



The reduction of the parallel thermal conductivity relative to the Spitzer value, caused by mag- 
netic mirrors, depends on the ratio Iq/X = 160 (Zo/lOkpc) (r/lO^K)"^ (^n/lQ-^cm-^) [Malyshkin & 
Kulsrud 2000]. The typical values are k\\/ks ~ 1/10 for the hot and low-density halo of a galaxy 
cluster, and K||/ks ~ 1 for the cluster central region. We see that for the halo, where Xq ~ 1 Mpc, 
Iq ^ 10 kpc, T ~ 5 X 10'' K and ra ~ 2 x 10~^cm~^, we have /3 ~ 20 and r ~ 10""^. Thus, according 
to figures 4, the heat conduction is unimportant there, while according to figures 3, the effective 
thermal conductivity is reduced by a factor of ~ 2000. At the same time, for the cluster central 
region we have Xq ~ 0.1 Mpc, Iq ~ 3 kpc, T ~ 2 x lO''' K, n ~ 10"^ cm"^, /? ~ 7 and r ~ 0.03, so the 
heat conduction is important there, and the effective thermal conductivity is reduced by a factor 
of ~ 200. This conclusion agrees with those of Tau (1995) and of Rosner and Tucker (1989). Thus, 
with certain parameters, such as these, thermal conduction can be important even with such a large 
reduction factor. Figures 4(a) and 4(b) give the direct temperature evolution without reference to 
reduction factors, and perhaps because of the sensitivity of the reduction factor to its definition, 
they are probably the more useful expressions of thermal evolution. 



Suginohara and Ostriker in their hydrodynamic simulations of galaxy cluster formation en- 



in our model are allowed to random walk in three dimensions, while Tribble considered a restricted one-dimensional 
random walk. 
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countered a "cooling catastrophe" , which appears as a steep non-reahstic rise in the density profile 
of the relaxed core of a galaxy cluster because of the fast cooling in the core (Suginohara & Ostriker 
1998). The heat conduction into the core was among solutions suggested by them. However, they 
believed that the conduction can be neglected if it is reduced by a factor of 30 or more. In this 
paper wc find that the thermal conduction is very important in the cluster central region, despite 
the fact that in our models it is reduced by a factor of ~ 200. Thus, the conduction should be 
included in hydrodynamic simulations of galaxy clusters, even when the reduction factor is large. 
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